Kinetic roughening, global quantities, and 
fluctuation-dissipation relations 

Yen-Liang Chou^ and Michel Pleimling^ 

"'^Department of Physics, Virginia Tech, Blacksburg, Virginia 24061-0435, USA 
E-maih ylchouSvt . edu, Michel . PleimlingSvt . edu 

Abstract. Growth processes and interface fluctuations can be studied through the 
properties of global quantities. We here discuss a global quantity that not only captures 
better the roughness of an interface than the widely studied surface width, but that is 
also directly conjugate to an experimentally accessible parameter, thereby allowing us 
to study in a consistent way the global response of the system to a global change of 
external conditions. Exploiting the full analyticity of the linear Edwards- Wilkinson and 
Mullins-Herring equations, we study in detail various two-time functions related to that 
quantity. This quantity fulfills the fluctuation-dissipation theorem when considering 
steady-state equilibrium fluctuations. 
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1. Introduction 

Due to its omnipresence in many fields in physics and engineering, kinetic roughening 
has attracted much attention over the years, see [H [21 El IH [5] for some earher reviews. 
In many instances the properties of growth processes can be understood by analysing 
rather simple Langevin equations. These Langevin equations, which can be either linear 
or non-linear, provide the theorists with a class of systems that they can study in 
a systematic way, using a variety of techniques. In addition, these Langevin equations 
have a wide range of applications, ranging from different growth processes to equilibrium 
step fluctuations at the surface of a crystal [SI [71 El [HI [I0| and kinetic smoothening of 
interfaces [TT] . 



The simplest Langevin equations studied in this context are the linear Edwards- 
Wilkinson (EW) [12] 



Here x is a rf- dimensional position vector located on the surface of the substrate, whereas 
for growth processes /i(x, t) is the height of a column at position x at time t|j 

These equations, which depend on the two parameters u (for the EW equation u is 
the surface tension or elastic constant) and D (the noise strength), are used to describe 
deposition with different relaxation mechanisms. For the EW equation the surface 
current responsible for the smoothening can be viewed to be due to the gravitational 
potential, whereas for the MH equation a surface current arises because of a chemical 
potential difference. A standard way to analyse kinetic roughening is to study the 
surface width, which typically displays the following three regimes. At the early stages 
the surface grows in an uncorrelated way during a regime that is sometimes called the 
random deposition regime. At a flrst crossover time ti, correlated growth sets in. This 
correlated regime continues until a second crossover time t2 at which the system enters 
the steady state regime. This time ^2 depends on the size of the system and shifts to 
larger values when increasing the linear extend of the substrate. For an inflnitely large 
substrate t2 diverges and the system never reaches the steady state. 

Quite some attention was paid recently to the ageing processes [E] that take place 
during the correlated growth regime [T5lll6l[I7lll8l[T9l|20l[2ll[221[23]. Most of the 
studies focused on local quantities as for example the height-height correlation function. 




(1) 



(2) 



(^(x,t)) = 



(r/(x, t)r]{x', t')) = D5\i^ - x')(5(t - t') . 



(3) 



X In the following we mostly use the language of growth processes, but due to the different physical 
situations described by the same Langevin equations our results have a broader range of applications. 



Kinetic roughening, global quantities, and fluctuation- dissipation relations 



3 



the response of the height to a local perturbation, the two-time roughness or the two- 
time incoherent scattering function [151 ESI EZl [13 [El [20]. In [21] we discussed the 
correlation function of the squared width and the response of the squared width to 
a global perturbation. Studying changes both in z/ and D, our work revealed that 
the global response of the surface width depends on how the system is perturbed. In 
addition, we showed that in the correlated regime the limit value of the corresponding 
fluctuation-dissipation ratio only yields the trivial value zero, due to the fact that the 
studied quantity, the square of the surface width, is not conjugate to any of the system 
parameters that are changed in our protocols. 

In this paper we propose to study a global quantity directly related to the effective 
Hamiltonian (m being an even number) 

H^ = yd'x{v-Hy (4) 

showing up in the Hamiltonian description that yields the corresponding stochastic 
equation of motion: 

Inserting @ into ([5|) yields the linear Langevin equation 

^^^ = -K^VrMx,t)+r^(x,t), (6) 

where for m = 2 we recover the EW equation, whereas m = 4 yields the MH equation. 

As we discuss in the following, our quantity, which can be computed exactly for 
linear Langevin equations, has many advantages. On the one hand, it better describes 
the roughness of a surface than the surface width itself. On the other hand, as this 
quantity is conjugate to u, a change in z/ yields a global response that allows us to study 
the corresponding fluctuation-dissipation ratio [211 [25]. For the case of an equilibrium 
steady state (as it is for example encountered for step fluctuations on crystal surfaces) 
the celebrated fluctuation-dissipation theorem is recovered, something that was not the 
case for the surface width [21] . 

The paper is organized in the following way. In Section 2 we introduce our quantity 
that we calculate exactly for the linear Langevin equations used in the context of non- 
equihbrum growth and related problems. Section 3 is devoted to the corresponding 
correlation and response functions. We thereby show that we recover the fluctuation- 
dissipation ratio for equilibrium steady states. Section 4 gives our conclusions. 

2. Global quantity conjugate to u 

Inspection of the effective Hamiltonian @ allows us to define both for the Edwards- 
Wilkinson and the Mullins-Herring equations the following time-dependent global 
quantity that is conjugate to u: 

GUt) = \jd'x{v^/'hy , (7) 
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For m = 2 the quantity Gm is of course readily identified with the total "kinetic energy." 
In the following we consider as substrate c?- dimensional lattices with linear extend L. 
For m = 2 our quantity can then be written as 

G2{t) = Wi: [Mx + a„ t) - /^(x, t)] ' (8) 



=1 

whereas for m = 4 we obtain 

d 



1 f 

= ^ E 1 E [H^ + ai, t) - 2/i(x, t) + /i(x - an, t) 

X U=l 



(9) 



The vector x now labels all points of the substrate lattice, whereas the vectors aj, 
i = 1, - ■ ■ ,d, indicate the primitive vectors on the substrate. For simplicity the lattice 
constant is assumed to be unity. Note that geometrically G2 is the total surface slope [§| 
and G4 is the sum of the curvatures. 

Usually the roughness of a growing interface of columns of hight /i(x) at substrate 
site X is expressed by the surface width 

w{t) = ^;^E(Mx,t)-Mt))' (10) 

where h(t) = jdY. h{yi, t) is the average height at time t. The quantity Gm can also 

be used to quantify the roughness of a surface [Jj]. In fact, it can even be argued that 
Gm captures the roughness much better than W. To see this, consider the two surfaces 
shown in Fig. [1] which have the same number of deposited particles. Whereas intuitively 
we would judge the surface (b) to be rougher than the surface (a), the surface width 
yields the same value for both cases. The value of Gm, on the other hand is larger for 
surface (b) than for surface (a), and this for both cases m = 2 and m = 4. 



(a) (b) 



Figure 1. Two surfaces with tlie same surface widtli but different values of Gm- 

For our purpose we need to derive the exact expressions for the average of Gm 
from the solution of the corresponding linear Langevin equations. Writing both the 

§ A generalisation of this quantity to the height difference between columns separated by a displacement 

r has been called the second order height difference correlation function l^j . 

II See j26j for an example where the surface slope G2 has been analysed in a numerical study 
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height /i(x, t) and the noise ?7(x, t) as a sum over reciprocal lattice vectors, we obtain 
the following expression 

(a^) = 2^^-^ y: {Kh-^)n^r'\ (n) 

where 

P(q)=X:(l-cosfe)). (12) 

1=1 

Here (■ ■ ■) indicates an average over the noise. The averaged quantity {Gm) is therefore 
a sum over two-point correlation functions weighted by the factor 'P(q)™'/^. Using the 
expression [21] 

{^(*.)"p(*2)> = ^«-""""*''"*"^;^ (^"""■■^'■""^ - 1) (13) 

for the two-point correlation function, with q = |q|, = |p|, and t< being the smaller 
of the times ti and ^2, we obtain 

= E 4 (i - ^^'"'"i n^r"- (14) 



The behaviour of (Gm) is therefore controlled by the length scale It = {2utY^"^, similar 
to the surface width [21]. Depending on the relation of It to the maximum and minimum 
values of g, qmax = iry/d and qmin = 2vr/L, different regimes can be discussed. 

1. For It < 1/qmax, the system is in the random deposition regime. An expansion 
in small It yields the expression 

(Gm)^2^Z}tEp(qr/2 



such that 



and 



^2'^Dtl^^j rfqP(qr/2 (15) 



(G2) ^ dL'^Dt, (16) 



{G4) ^d{2d + l)L'^Dt. (17) 

In this regime (Gm) varies linearly in time and shows the same time dependence as the 
squared surface width [21]. 

2. When It > l/qmin, the system is in the saturation regime and (Gm) reaches its 
maximal value, 

(Gm)=G,m-2^:^E4p(q)-/2 

D /L^ 



m-4 D f L\ 

2— _ _ I{m,d), 

V VTT/ 



with I(m,d) = Jq 'P{q)"^^'^ /q'^dq. The numerically evaluated values of I{m,d) for 
various d and m are shown in Table ([T]). 
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d=l 


d=2 


d=3 


d=4 


m=2 


1.21532 


3.43997 


10.29687 


31.53744 


m=4 


0.49733 


1.26083 


3.55376 


10.51476 



Table 1. The numerically evaluated values oi I{m,d) = JJ^ 7'(q)™/^/q™(iq 




t t 



Figure 2. (Gm) for (a) the Edwards- Wilkinson case m — 2 and (b) the Mullins- 
Herring case to = 4. Systems of various lengths and various dimensionality of the 
substrate are shown. The insets display the power-law approach to the stationary 
value Gs,m- The symbols are obtained by numerically evaluating the exact expression 
Eq. (HH). The dashed and solid lines are given by the asymptotic expressions (|16|17p 
and (|20|) . The system parameters are ly — 0.001 and D = 1. 



3. To calculate the asymptotic behavior in the correlated regime, where 1/qmax < 
It < ^/qmin, we rewrite Eq. ( IT^ in the form 

and take advantage of the hyperspherical symmetry of 'P(q) for g < 1: P(q) ~ 
In fact, since ( fT9l) involves the factor e"''™'*", we can neglect contributions with g > 1, 
provided that It ^ 1. In that case we can simply replace the sum by an integral, yielding 
the result 

^Gs,m-^(^) nd^^{2vt)-^ (20) 
4z/ \27r/ m 

for the correlated regime, where VLd is the surface of the d— dimensional unit sphere. 

Figure |2] shows (Gm), obtained by numerically evaluating Eq. ( IT4l) . for various 
values of m and d and various system sizes. Also shown as lines are the asymptotic 
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expressions (ITBI) and (IT7I) (dashed lines) and (120]) (full lines). The plots and the 
asymptotic expressions show that (Gm) grows linearly with time t when the system 
is in the random deposition regime, the pre-factor being independent of the value of i'. 
After the system passes the first crossover point, a pre-factor that depends on z/ shows 
up. In the correlated regime {Gm) approaches the steady-state value Gs,m algebraically, 
with an exponent — . 

We should point out that the plots of {Gm) versus time, shown in figure [2l do 
not easily reveal the crossovers between the different growth regimes, in contrast to the 
surface width where the different regimes are readily identified by mere inspection. This 
might somehow restrict the usefulness of (Gm) , especially in situations where the steady 
state value is not readily known. Still, as we discuss in the following, the quantity Gm has 
many other advantages that makes it appealing both for theoretical and experimental 
studies. 

3. Correlation function, response function, and fluctuation-dissipation ratio 

Two-time quantities are extremely useful when studying relaxation phenomena, as they 
allow to capture most of the processes that underly the properties of a system far 
from its steady state. Typical quantities are the correlation and response functions, 
as well as combinations of these as for example the fluctuation-dissipation ratio. In 
the context of non-equilibrium growth processes and interface fluctuations both local 
p3| [T6| \T7\ [T8| [T9| [20| [TT| and global [21] quantities have been studied to some extent. 
The former include the height-height correlation function, whereas the change of the 
surface width to a global perturbation is an example of the latter. In the following we 
are deriving the corresponding exact expressions that involve the global quantity Gm- 

We first consider the two-time correlation function of Gm- Starting from Eq. ( ITTi) 
we have 

{Gm{t)Gm{s)) = 2^-'L'' Vic^r'^VivT'^ 

{{h^{t)h_^{1:)){h^{s)h_^{s)) 
+ {h^{t)h^{s)){h_^{t)h_^{s)) 

+ {h^{t)h^^{s)){h^^{t)h^{s))}. (21) 

Here we note that the first term yields {Gmit)){Gm{s)). Since V{q) is an even function, 
we have 

P(q)™/^P(-q)"/2 = P(q)™/2p(q)'"/2 = P(q)"^. (22) 

Using the Eqs. ( fT3|) and ( !22l) we obtain the connected two-point correlation function of 
the quantity Gm '- 

GG{t,s) = {Gmit)GUs)) - {Gmit)){Gmis)) 

= 2— 3:^^ J2 -^e'^'"^"' {e''"^"' + e-^'^"' - 2) P(q)'". (23) 
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We can also calculate the evolution of the average of Gm subjected to a perturbation 
when we suddenly change u during the growth process. This is a natural thing to do as 
in the effective Hamiltonian (jl]) u and Gm are conjugate quantities. Note that this is 
the response of a global quantity to a global perturbation. 

In this case, the evolution of the average of Gm can be written as 

{Gm),-.At, s) = y: {h^r''^h^'^^'^)v{c^r'\ (24) 

where the notation "/i — )■ z^" indicates that the change from to at the waiting time 
s. 

As discussed in [2T1 [20], the solution of the Langevin equation at times t > s 
becomes under that change 

/,(^^->-)(t) = e-'^'i'"(*-^)/iq,^(s) + r dre-'^'''"(*--)r/q(r), (25) 

J s 

where 

h^As)= rrfre-'^'''"(^--V(r) (26) 



is the solution of a surface that evolves until time s at the value fi when starting from 
a fiat initial state. 

Plugging Eq. ( I25l) into Eq. ( 12^ and comparing with Eq. ( HM , we straightforwardly 
obtain the response function 

{Gm) {t,s)-{Gm)u{t) 



XG{t,s) 



e 



:i-e-^'"-)jp(q)f, (27) 



1/ 

where e = z/ — /i. 

With these exact expressions for the correlation and response functions, we can now 
discuss the properties of the two-time quantities in the different regimes. In addition 
we can also study the behavior of composite quantities, as for example the fluctuation- 
dissipation ratio: 

Xit,s)^^^^/^^4^. (28) 
as OS 

The asymptotic steady-state behavior (i.e. the large us limit) is easy to obtain for 
both quantities. Considering only the contribution from terms with minimum q, we 
have 

Cc(t,s)~^e-2C«"'^(*-^) (29) 

and 

XG(t,.)^^e-2««^'^(*-^), (30) 
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where we have apphed the approximation 

nqmin)/qLn ^ 1/2 (31) 

to both expressions and, in addition, the hmit e = z/ — /i— T-Oto the response function. 
For the fluctuation-dissipation ratio (!28|) we obtain in the steady state, with t — > oo, 

X = ^. (32) 

Assuming the validity of the Einstein relation D = 2T, which makes the steady state to 
be an equilibrium steady state, we recover the fluctuation-dissipation theorem X = 1/T, 
as expected. 

The most interesting case is the case where both the waiting and observation times 
are in the correlated regime (to leading order identical results are obtained when the 
waiting time is still in the initial RD regime). This corresponds to l/qmax < h,h < 
1/qmin, with Is < It and ^ 1. In this limit, we replace in both Equations fl25]) and fl27|) 
the sums by integrals and treat the integrands as hyperspherical symmetric functions 
(this is the same method we used to derive the asymptotic equation for {Gm) in the 
correlated regime). We thereby obtain the following power-law decay functions 



CG{t,s)^Dh^-j n,^^-^{2ut)-'-^, (33) 



XG{t, s) ^ Ds' ( + [2uit - s)]-'-^ (34) 



and 



.2nJ 2m 

where we replaced e^''"'^^ + e"^'^'"'''' — 2 by (2g''"z/s)^ for the correlation function and 
(1 — e~'^''^^^)/ ^ ~ (1 ^ e"^^""^*)/// by 2{v — /i)s^g^™' for the response function before 
integrating. These replacements correspond to retaining only the leading terms in the 
Taylor expansions. 

In the literature on physical ageing it is convention to write for a system undergoing 
simple ageing the two-time correlation and integrated response functions in the form 



Ccit, s) = s-'fcit, s) ; xcit, s) = s-'^f^it, s) , (35) 

where fc{y) and f^iy) are scaling functions which decay algebraically for large 
arguments: 

/c(y) ~ y-'^/^ ; fxiy) - y-'''^^ (36) 

where Ac and are the autocorrelation and autoresponse exponents, whereas z is 
the dynamical exponent (for our linear Langevin equations we have that z = m). 
Recasting Equations (!33l) and in these ageing forms, we immediately obtain that 
a = b = d/m, whereas the autocorrelation and autoresponse exponents are given by 
= ~ 2 + d/m. This ageing scaling is illustrated in Figure [3] for the one-dimensional 
EW system. 
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Figure 3. Ageing scaling of (a) the correlation and (b) the response function of G2 
in the one-dimensional EW system. The parameters for the calculations are L = 2^^, 
f = 0.001, and D = 1. 



Let us close this Section with a more careful discussion of the fluctuation-dissipation 
ratio (128|) . As mentioned in the introduction, in a finite system the correlated regime 
goes over into the stationary state at a system size dependent time ^2- In the infinite 
system ^2 diverges and the correlated regime prevails for all times. For the limit value 

Xoo = lim lim X(t,s) (37) 

s — >oo t — >oo 

we therefore obtain the values X^o = 1/-D = 1/2T for the infinite system and 
Xoo = 2/D = 1/T for the finite system, where we introduced temperature via the 
Einstein relation. The finite system ending up in the steady state for finite times, we 
recover the fluctuation-dissipation ratio. If the system remains in the correlated region, 
the effective temperature is twice that of the heat-bath. The crossover between these 
two regimes can be visualized for finite systems by plotting X{s) for t ^ s as done in 
Fig. m For that figure we plot the value of X{s) = X{s + 10^, s) as a function of s, 
which yields the value X{s) = 1/D for s ^ t2 and the value X{s) = 2/D for s ^ ^2- 
The crossover times for m = 2 and m = 4 are indicated in Fig. H] by the vertical lines. 
Fig. Ogives a more comprehensive view of the behaviour of X{t, s) as a function of both 
s and t for the one-dimensional EW equation. Two plateaus can be distinguished in the 
contour plot: one for the steady state (i.e. the regime where t > t2 and s > t2) where 
X = 2/ D and one away from stationarity, with t^ l/2v and s < t2, where X = 1/D. 
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Figure 4. The fluctuation-dissipation ratio for m = 2 and to = 4 in various dimensions 
as a function of the waiting time s. Here X{s) — X{s + 10®, s). The hnear extension 
for afl systems is L — 2^. The full (dashed) line indicates the crossover time t2 for the 
EW (MH) equation. Note that the value of X{s) is independent of the dimensionality 
of the substrate. 



4. Conclusion 

In the past the study of global quantities in systems relaxing towards a steady state 
has proven very fruitful in a large variety of systems (see the corresponding discussion 
in [Hj). In [2T] we did a first attempt at using global quantities in the context of 
correlated growth and interface fluctuations, choosing the surface width as our global 
quantity. However, the surface width is a complicated quantity that has the notable 
drawback that the conjugate system parameter is unknown. Consequently, it is not 
possible to form a meaningful fluctuation-dissipation ratio using that quantity. 

In this paper we are proposing a different global quantity for the study of kinetic 
roughening and related interface problems. This quantity is proportional to the effective 
Hamiltonian used in the Langevin description and is conjugate to a system parameter 
that can be changed in experiments [201 [11]. f^ict, Gm seems better suited to capture 
the roughness of a surface than the surface width itself, as illustrated in figure [H 

Focusing on linear Langevin equations we derive exact expressions for Gm as well as 
for the corresponding correlation and response functions. This allows us to discuss also 
more complicated quantities as for example the fluctuation-dissipation ratio. In fact, 
we recover for the quantities derived from Gm the fluctuation-dissipation theorem for 
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1 Ti2 Tx3 4 T25 6 7 

log(t) 

Figure 5. Contour plot of X{t,s) for the one-dimensional EW equation. The 
parameters used here are L — 2^ , v = 0.001, and /i = 0.99;^. Also indicated are the 
logarithms of various relevant time scales: n = logii, T2 = logt2, and Tx — log(l/2i^). 
Here and in the following figures, we illustrate our results for the one-dimensional 
EW equation, but, as the exact results reveal, similar results are obtained for the MH 
equation as well as for substrates of higher dimensionality. 

equilibrium steady states, whereas in the correlated regime we can assign an effective 
temperature to our system. 

All calculations presented in this paper have been done in the context of linear 
Langevin equations. However, most growth processes are governed by non-linearities. It 
is therefore important to clarify to what extend our results obtained for linear Langevin 
equations remain valid when considering non-linear stochastic equations as for example 
the Kardar-Parisi-Zhang equations [27] . We intend to address this and other questions 
in the future. 
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